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Abstract. Following the method of Hoenselaers and Perjes we present a 
new corrected and dimensionally consistent set of multipole gravitational and 
electromagnetic moments for stationary axisymmetric spacetimes. Furthermore, 
we use our results to compute the multipole moments, both gravitational and 
electromagnetic, of a Kerr-Newman black hole. 



1. Introduction 

After Ernst's first paper on reformulating Einstein's equations for stationary and 
axisymmetric spacetimes (i) without PP, and (ii) with electromagnetic field [2], 
Geroch and Hansen ^ E] defined the multipole moments of static and stationary 
spacetimes respectively, by measuring the deviation of the geometry from flatness in 
the neighborhood of a point at infinity. Their definition of the moments is equivalent 
to the one introduced by Thorne in Later, Simon and Beig [HI reconstructed the 
moments from the series expansion of field variables, and Simon [7] generalized them so 
as to include electrovacuum spacetimes. Especially the axisymmetric spacetimes that 
Ernst studied are interesting physically, and can be fully described by a number of scalar 
moments, in the same way that newtonian gravitational potentials are characterized by 
their mass multipole moments. 

In 1990 Hoenselaers and Perjes [S] generalized the work of Fodor et al ^ on vacuum 
gravitational moments to cases where spacetime is endowed with an electromagnetic 
field that retains the axisymmetry of the metric. Although the generalization method 
they used was right, an incorrect intermediate formula led to the contamination of 
the final expressions for the moments with errors. In our paper we have tried to (i) 
determine the root of these errors, (ii) give a few extra intermediate relations that are 
used in the procedure of computing the moments, and (iii) correct the final expressions 
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that relate the gravitational and electromagnetic multipole moments of spacetime with 
the power series coefficients of the gravitational and electromagnetic potentials. The 
final expressions for the gravitational and electromagnetic moments may be important 
for people that are looking for approximate expressions for the metric of stationary 
and axisymmetric sources that depend on a few physical parameters of the source. We 
actually needed these expressions when we tried to extend [TU] Ryan's study (TT] on ways 
to determine the moments of a compact object from gravitational wave data analysis, to 
axisymmetric spacetimes with electromagnetic fields. The dimensional inconsistency of 
the moments of jH] , that we were to use, drew our attention and we set out to thoroughly 
check the whole procedure that Hoenselaers and Perjes followed in their paper. 

The rest of the paper is organized as follows. In Section|21we briefly repeat the basic 
steps of the method used, giving a few extra intermediate expressions that may be useful 
to the reader. Although most intermediate steps can be found in a number of previous 
articles, the plethora of different conventions and notations used by different authors, 
urged us to systematically present the ones we used to calculate the moments. In Section 
Elwe present the ffist five gravitational and electromagnetic moments as functions of the 
power series coefficients of the gravitational and electromagnetic potentials. A number 
of comments on the terms that show up at each moment have been added at the end. 
In Section 0] we compute once again the ffist five power series coefficients of the mass 
and electromagnetic potentials of a pole-dipole source as in |Hj. Finally, in Sectional 
we use the whole technique to compute the moments of the Kerr- Newman metric, from 
the knowledge of the metric along its axis of symmetry. The gravitational moments are 
exactly like in Kerr metric, while the electromagnetic moments are simply e/m times 
the corresponding gravitational moments. 



2. Ernst functions and moments 



The metric of a stationary axisymmetric spacetime could be described by Papapetrou's 
line element 



ds^ = -F{dt - ujd(Pf + F-^ [e^^dp^ + dz^) + p'^dcj)'^] , (1) 

where -F, cu, and 7 are functions of p and z. Assuming that spacetime is asymptotically 
flat, 7 — > 0, F ^ 1, and u; ^ at infinity. As Ernst has shown, any axially 
symmetric solution of the coupled Einstein - Maxwell equations can be described by 
two complex functions £ and $ which are related with the metric functions F and u 
and the components of the electromagnetic potential At, A^ (cf., Eqs. (8,9,11,14) of 
The third metric function 7 can be deduced by using the two remaining Einstein 
equations. 

Now, instead of £ and $ we use the complex Ernst potentials ^, g, which are the 
analogues of newtonian gravitational potential and Coulomb potential respectively, 

^ = ^, $ = ^. (2) 

1+e 1+e 
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Note that the former relation is the opposite of Ernst (cf., [1112]), and thus our ^ is the 
inverse of Ernst's C,e- It is easy to show that the Einstein equation for S and $ [8; take 
the following form for ^ and q 

- qq* - = 2(rVe - q*Vq) ■ (3) 

- qq* - l)V2g = 2(rve - q*Vq) ■ Vq. (4) 
By * we denote the complex conjugate. The Ricci tensor of the three-metric is given by 

m 

(a* - qq* - lYR^J = 23? (V,^V,r - V,gV,g* + s,s*) , (5) 

where 

S^ = ^V^q - qV^^. (6) 

Vj denotes the components of the gradient, and 3fJ() denotes the real part. 

In order to measure the moments of an asymptotically flat spacetime, according to 
Geroch-Hansen procedure [Sill], we map the initial three- metric to a conformal one 

hij > h^j h^j. (T) 

The conformal factor Q should satisfy the following conditions: fi|A = DiQ\\ = and 
DiDjfl\/^ = 2hij\\, where A is the point added to the initial manifold that represents 
infinity. Q transforms ^ and q potentials to 

I = q = Q-y\ (8) 

By a coordinate transformation 

P 9 I 9 ' 



= (9) 

we bring infinity at the origin of the axes (p, z) = (0, 0). Then, by choosing the conformal 
factor to be 

= f2 = p2 + (10) 
the conformal metric in the new coordinates takes the following form 
/ e^^ \ 

hij = e^^ , (11) 
V J 

which is flat at f = 0, since 7|f=o = 0. The Ricci tensor of this conformally transformed 
three-geometry in barred coordinates is 

(r'rl - r'q*q - l fR^J = ^^mD^C - D,~qDfq* + 5,5*), (12) 
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as one could verify by a tedious but straightforward calculation, with 

dp dz 
9 d ^ 

D2 = PT- + Z— + 1 

op oz 

~s, =f{lD4-qDi). (13) 

The DiS are the Vj's of jB^ with their indices interchanged. Notice that this has exactly 
the same form as Eq. (jS)) if we replace the free (not differentiated) ^'s and g's with and 
fg, respectively, while the differentiated ^'s and g's simply with ^ and g, respectively. 
In Eq. (6) of |H] which was used to evaluate the Sj's, the free ^ and q were replaced by ^ 
and g, instead of and fg, to get the Sj's. This is exactly what caused the dimensional 
inconsistency in the final expressions for the moments in jSj. The field equations ()3|4j] 
keep exactly their form in the conformal metric if expressed in the barred coordinates, 
by using the same replacements. 

The tensorial multipole moments of the electrovacuum spacetime are computed by 
the recursive relation (see [Oj) 
p(o) = I 

P^' = 

P&l,. = C[V.„,,iti - ln{2n - l)i?../t^:ij, (14) 



for the geometry, and 

Q(°) = g 
= I 
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for the electromagnetic field. The symbol V is used to denote the covariant derivative in 
the conformal space and it should not be confused with the same symbol used in jH] . The 
operator C denotes the operation "symmetrize over all free indices and take the trace-free 
part". Note that all tensors should be evaluated at A (at infinity). Due to axisymmetry, 
the components of these tensorial moments are multiples of the corresponding scalar 
moments which are the projection of the tensorial moments on the axis of symmetry. 
Following Fodor et al jH], the scalar moments are defined as 

1 

• 9 

A 

(16) 

A 

Note that these moments have the opposite sign from the ones defined by Hansen 
The gravitational scalar moments are then simply given by 

(2^-1)!!^°"^' ^^^^ 



n ~ 1 -'2.. .2 

n! 

Qn — r Q2...2 
nl 



On the Multipole Moments of Axisymmetric Electrovacuum Spacetimes 5 

where S^^ are computed by the recursive relations (23) of jSj- The electromagnetic 
scalar moments Qn are computed from exactly the same recursive formulae, if we replace 
the initial term S^'' with q instead of In order to use these recursive formulae one 
needs also the derivatives of 7 which could be expressed as linear combinations of the 
components of the Ricci tensor which is given in Eq. (fT^ 

71 = 7,p = ^p{Rpp- R-zz), 

72 = l,z = pRpz- (18) 

We end up this section by writing both ^ and q as power series of p,z: 

00 

i,j=0 



00 



q=}_^bi,p^z=, (19) 

i,j=0 

where aij and bij vanish when i is odd (this reflects the analyticity of both potentials 
at the axis of symmetry). Because of Einstein's field equations (see Eqs. fl3l4|) and the 
discussion for the tilded version of them which follows after Eq. p3j) ). the coefficients 
in the above power series are algebraically interrelated: 

(r + 2f ar+2,s = - (s + 2)(s + l)ar,s+2 + 

+ ^^kl^*mn - bklb*mn) X 

k,l,m,n,p,g 

X [apg{p^ + - Ap - hg - 2pk - 2gl - 2) + 
+ ap+2,g-2{p + 2){p + 2-2k) + 

+ ap^2,g+2ig + 2){g + l~2l)], (20) 



and 



(r + 2f br+2,s = -{S + 2){S+ l)br,s+2 + 



V 



+ ('^kl<Knn ~ bklbl 

k,l,m,n,p,g 

X [bpg{P^ + g^-Ap-5g- 2pk - 2gl - 2) + 
-t- «p+2,g-2 {p + 2){p + 2-2k) + 

+ b,^2,g+2{g + 2){g + 1-21)1 (21) 

where m = r — k — p,0<k<r,0<p<r — k, with k and p even, and n = s — I — g, 
< / < s + 1 , and —l<g<s — l. These recursive relations could built the whole 

power series of ^ and q from their value on the axis of symmetry 

00 00 

ap = 0) = 5^m,^\ g(p = 0) = (22) 

i=0 i=0 

Thus from ^(p = 0) and q{p = 0), we can read the multipole moments of spacetime. In 
the next section we will write the expressions that connect the multipole moments of 
both types with the m^'s and gj's. 
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3. The multipole moments of stationary axisymmetric electrovacuum 
spacetimes 

Following the procedure that is outlined in Section |21 we compute the first five 
gravitational and electromagnetic multipole moments in terms of aij and bij. Then 
the algebraic relations between aij and bij are used in order to express the moments in 
terms of mj = aoi and qi = b^i. If we define the following useful quantities 

Mij = miirtj — nii^inij^i, 

Qij = QiQj — Qi-lQj+li 

Sij = miQj - mj_igj+i, 

Hij = qiUij - gi_imj+i (23) 

the gravitational moments as functions of the power series coefficients of ^ along the 
symmetry axis are given by 





= mo, 


Pi 


= fUi, 




= m2, 




= fns, 


Pa 


= m4 


P5 


= fn^ 



1 13 

1114 
-m*M3o - — m*M2o + -QoSso + ^^o-^si 



- ^QlSu - ^q;S,o, (24) 
while the electromagnetic moments as functions of the power series coefficients of q along 



axis 


are g 


Qo = 


go, 


Qi = 




Q2 = 




Q3 = 


13, 


Qa = 


qA + 


Q5 = 


^5 + 



1 3 

1114 
g^oQso + ^gt<520 - ^rnlH-iQ - —mlH2i + 

1 1 

+ i^rn\Hu + —mlHiQ. (25) 
One could make the following comments on the expressions above: 

(i) While the quadratic quantities Mij, Qij are independent quantities, Sij, Hij are not 
completely independent of each other. The following obvious relations between 
them hold 

Sij = —Hj^i (26) 
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Thus, although the introduction of four quantities looks like a redundancy this gives 
a symmetric form in the two sets of moments. 

(ii) The expressions for the moments coincide with the expressions of [H], if there is no 
electromagnetic field. 

(iii) The different expressions we have in -P4, (54, -P5, moments is due to incorrect 
replacement of q,^ with g, ^ in Si in Eq. 0. As was noted in Section 121 the new 
expressions for the moments are now dimensionally right (see note (v) below). 

(iv) From the form of the expressions we notice that there is a symmetry between 
gravitational and electromagnetic moments. The electromagnetic moments arise 
from the corresponding gravitational moments if we interchange qi with rrij, replace 
Qij with Mij and Sij with Hij, and reverse the signs of all but the single rrij, qi 
terms. 

(v) Since ^ should be dimensionless, and qi should have dimensions [M]*"*"^, and 
hence all terms in each moment have the same (sum of indices) + (number of 
indices). This simple observation gives us the capability of predicting what kind 
of terms we expect at each moment. This dimensional analysis makes it clear that 
the expressions (27,28) of [S] contain incorrect terms. 

4. A spinning mass with charge and magnetic dipole 

For completeness, we rewrite the results of the pole-dipole example used by Hoenselaers 
and Perjes, based on the correct expressions for the moments. Thus, we also assume 
that the moments of the source are 



Pq = m, Qo = e. 

Pi = iam, Qi = i/ie, 

Pn = 0, and Qn = 0, for n > 2. 



(27) 



Then, if we invert the expressions of Eqs. ()24p25|l we get 



mo 



m 







7713 









and 



go = e 



qi = ifie 



q2 = 
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?3 = 

Q4 — -lJ'e{am'^ — jie^) + —eam^{a — jj) 



Q5 = 



/ 2 2 2 2\ 

— //e(a m - II e) 



(29) 



Although 777-2 — 0.2 — 'm'3 = ?3 = 0, all higher order and Qi do not vanish. Especially, 
if we set — a, 

mo = m 
mi = iam 
m2 = 
ms = 



2/2 2\ 

TTT'A — -^jO- m[m — e ) 



and 



7715 — — t 



go = e 

qi — iae 

52 = 

?3 = 
1 



1 

21 



a^m{m'^ — e^) 



(30) 



-^2/2 2\ 

?4 = e(7?7 — e ) 



55 = 



21 



a e[m — e ) 



(31) 



which are exactly what the formulae of Hoenselaers and Perjes yield. If on the other 
hand we set = 777^, then 

777.0 = m 

mi — iam 

7772 = 
7713 = 



1 3 

-a H 

7 70 

1 



7714 = "^^(^fl + -^jj,){a — jj) 



and 



7715 = - 

qi = ifie 



21 



am {a — n ) 



(32) 
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q4 = e\-fx + —a){a- fx) 



= - « 



—fie (a - /i J 



(33) 



which do not coincide with the corresponding expressions of Hoenselaers and Perjes. Of 
course if = and /i = a, then all qi and except the first pair of them vanish. 

5. The Kerr-Newman metric 

In this section we follow the technique presented here to compute the gravitational and 
electromagnetic moments of a Kerr-Newman black hole. Ernst constructs this metric 
in prolate spheroidal coordinates {x,y), by setting 

x + iay e 

iE = , qs = — • (34) 

m m 

The ^ we have used in our paper (and in the majority of relevant papers) is the inverse 
of Ernst's (compare Eq. (0) of the present paper with Eq. (10) of |.J|), and q = qE^- 
Additionally, we will use the opposite sign for a from the one used by Ernst, since his 
convention leads to a spinning black hole with its spin along the —z direction for positive 
a (see footnote 24 of Israel ^|). Thus, 

i = — , q = — • 35 

X — lay X — lay 

At the axis of the black hole (p = p = 0), the prolate spheroidal coordinates are y = 1, 
and X = z. Therefore, the Ernst potentials along the axis take the following form 

e(p = 0) = - = ^^, g(p = 0) = -e(p = 0), (36) 

z — la 1 — laz m 

and finally the conformally transformed potentials are 

oo 

e(p = 0) = - — = mY,{ 



, [lazY 

1 — iaz 

1=0 



oo 



q{p = 0) = — = eJ2i^azy. (37) 



laz 

1=0 



Now it is straightforward to read the multipole moments of such a source. It is easy to 
see that in this case, as in the case of Kerr metric, Mij = Qij = Sij = Hij = 0, so the 
mi's and g^'s coincide with the actual moments. 

Pi = m{ia)\ 

Qi = e{iay. (38) 

This set of moments fully defines the gravitational and electromagnetic field of a Kerr- 
Newman black hole. 
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